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1 Itroduction 



Powerful bipolar outflows appear to be common features of the wide class 
of astrophysical objects including protostars, young main-sequence stars, 
galactic X-ray sources, Active Galactic Nuclei (AGN), and quasars [1]. In 
the majority of these cases there exists more or less convincing observa- 
tional evidence about the presence of large-scale magnetic fields in these 
objects. Ordered magnetic fields should have the crucial role in giving rise 
to such outflows or jets. Namely, they may be responsible for the collima- 
tion of the jets and/or acceleration of matter up to relativistic velocities in 
them [2]. 

High compactness and huge energy output, which characterizes some 
(maybe, the most interesting) concrete kinds of objects featuring bipolar 
(or, unipolar) outflows (such as AGNs and quasars), may be explained if 
one assumes that they are produced in result of accretion onto the rotat- 
ing supermassive black hole. It seems evident that a consideration of the 
innermost part of such outflows, in the close neighbourhood of the black 
hole must be performed in the framework of general-relativistic magneto- 
hydrodynamics (MHD). Such treatment enables us to take into account 
properly the influence of strong gravitational and electromagnetic fields 
onto the structure of the jet. Since 3 + 1 formulation of black hole elec- 
trodynamics [3-5] is the most convenient mathematical apparatus for such 
purposes, in the present study we shall establish our consideration on it. 

It must be noted from the very beginning that the method, which is 
developed in this paper is the generalization of the advanced methodology 
used in the number of recent references for non-relativistic [6-9], and for 
special-relativistic [10-11] magnetically driven jets and winds. By means of 
the method we derive the set of equations for general-relativistic jets from 
the equations of the ideal 3+1 MHD. Under certain simplifying assump- 
tions we obtain the simple, representative solution of these equations. The 
features of the solution are discussed and compared with the one, which 
has been found in [8]. 
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2 Governing equations 



In the forthcoming analysis we shall use the following notations: (a) greek 
indices will range over t, r, 9, (f) and represent space-time coordinates, com- 
ponents, etc.; (b) Latin indices will range over r, 6, (f) and represent coordi- 
nates in three-dimensional "absolute" space. We will approve the spacelike 
signature convention (— + + +) and units in which the gravitational 
constant and the speed of light are equal to one. 

The rotation of a central object (i.e. a rapidly rotating Kerr black 
hole) introduces off-diagonal terms g t ^ in the metric so that the space-time 
generated by the rotating body is represented by the metric: 

ds 2 = g tt dt 2 + 2g t(f) dtd(f) + g^dcj) 2 + g rr dr 2 + geedO 2 : , (1) 

with the metric coefficients independent of t and (f). 
In 3+1 notations (1) may be rewritten as [3-5]: 

ds 2 = -a 2 dt 2 + j lk {dx l + (3 l di)(dx k + (3 k dt), (2) 

where a is the so-called lapse function defined as: 

2 

a , (3) 

9<M> 

7ifc is the three-dimensional "absolute" space metric tensor (with diagonal 
nonzero components ^^=7^) and (3 is the spatial (three-dimensional) vector 
with components 

/?Wo, 0, ^\ (3i = j ik (3 k . (4) 
V 9<M>) 

Note that the Kerr metric is the subclass of the general metric (1). In 
particular, 

(■> 2Mr \ Sr^ 

gtt = -[l-—y (5) 
2aMrsin 2 9 

9t4> = ^ , (6) 

Asin 2 9 

9H = — £ — j ( 7 ) 
S 

9rr — X' \°) 
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gee = S; (9) 

where a=J/M is the specific angular momentum of the black hole per its 
unit mass, and 

E=r 2 + a 2 cos 2 6, (10) 

A=r 2 -2Mr + a 2 , (11) 

A=(r 2 + a 2 ) 2 -a 2 Asin 2 9. (12) 

In the present study we consider jets as being axisymmetric {3$ = 0) 
and stationary (d t = 0). We, also, neglect dissipative effects. The ba- 
sic equations are those of the general-relativistic MHD written in 3 + 1 
formalism: 

V[Tmn(av - (3)} = 0, (13) 

Vx(aB) = 47raJ - C$E, (14) 

E + vxB = 0, (15) 

Vx(aE) = C^B, (16) 

VB = 0, (17) 

rip 1 

f- = -~ 2 V(a 2 S) - ^T l \ (18) 

r /c 1 

— = e^+H-5--V(aT). (19) 
dr a 

In these equations we use the following notations: 

r = [1 - (v-v)}- 1 ' 2 , (20) 

iL=I [ft _(£ V )], (21) 

^-±V(a), (22) 
cm 

Hik=—/3k-.i, (23) 
en 

o"zfc= — + ^*ife); (24) 

and all vector and tensor quantities are defined in the three-dimensional 
"absolute" space with the metric 7^. 
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Note also that in these equations n is the proper baryon number density, 
while E, B, and J are the vectors of the electric field, the magnetic field 
and the current density, respectively. The "Lie derivative" of a vector A 

— * 

along the vector (3 is defined in the following way: 



C ? A=(A, V)/J - 03, V) A, 



(25) 



and v is 3-velocity, related to the spatial components of matter 4- velocity 
w a =(V, u l ) via the expressions: 



i 1 ( Ul oi 
In (18) e is the total energy density, defined as [5]: 

£=£p + £f, 

£ p =(mn + Pv 2 )T 2 , 

5 is the total momentum density, defined as: 

S —Sp -\- Sj, 
S p =(mn + P)T 2 v, 



and Tik is defined as 



T f =— 

1 ik— ■ 



4tt 



T^(mn + P)r 2 W + P^, 

+ BiB k ) + ^(£ 2 + 5 2 )^ fc 



to: 



Note that the roots of the determinants of \\g a p\\ and 

^/g=(det{g a! 3)) 1/2 = EsinO, 



(26) 



hk 



(27) 
(28) 
(29) 

(30) 
(31) 
(32) 

(33) 
(34) 

(35) 

are equal 

(36) 
(37) 
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3 Main consideration 



Let 9(r) denote the jet outer boundary and 



2 

0) ' 



(38) 



be the dimensionless angular variable. From now on, we shall consider 
"narrow jet" i.e., 9(r) is assumed to be as small that sin9~9 and cos9~l. 
Surely, the same is true for the angular coordinate 9 < 9{r) inside the jet. 
The assumption noticeably simplifies the angular dependence of various 
quantities appearing in the theory. For example, from (37) we get: 



r + a 



2\3/2 



and owing to the same assumption: 



A 1 / 2 



A 



-0, 



(39) 



1/2 



j 2 + a 2 J 
2aMr 



(40) 



(r 2 + a 2 ) 2 ' ^ 

Let us introduce, now, the magnetic flux function ip(r : 9) [8], defined as 
to have: 

B r = -^(e), (42) 

B e = ~d r m, (43) 

then it is easy to check out that the "no monopoles" condition (17) is 
satisfied automatically, and also 



de/dr 



or 



de/d9 

Inserting (44) into the (17) and resolving it for B r we get: 



(44) 



(r 2 + a 2 )<? 2 (r) ' 



(45) 



where B r is the "physical" radial component of the magnetic field vector, 
while the $ is some constant of integration. 
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Equations (15) and (16) may be combined to give the so called "nduction 
equation" : 

rot[(av- P)x§] = 0. (46) 

Poloidal components of this equation lead to the simple relation between 
poloidal components of v and B: 

v r B° = v 6 B\ (47) 

which, together with (44) leads to 

v o = jgMj. (48) 
or 

The toroidal component of (46), after taking into account of the (17) 
may be written as 

d r [a^v r B*} + d e [ay/TV 9 B+\ = + (49) 

In this paper we assume that Q = Q(r). It allows us to neglect the last 
term on the right in (49). Remained equation may be solved separately as 
inside [at 6 < 9{r)] as outside the jet [when 9 > 9{r)]. In such a way we 
get: 

B"{r)= 9<l (50a) 

BHr)= f ? {n ~ n i v 9> 9. (506) 
v ' (r 2 + a 2 )v r 9 2 {r) v ' 

Let us consider, now, the momentum conservation equation (18). For 
stationary case it may be rewritten as: 

--0-V)S = eg + U S - -V(af ). (51) 

a a 

For the 0-component of this equation taking into account that = 
and 

—P k S ( j ); k = H ( f,kS k , (52) 

we get simple equation: 

K/77j),* = 0. (53) 
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For the poloidal components of the same equation, taking into account 
that 

-P k S iik = -H 4>i S+ i = r,6, (54) 
a 

we derive the following equation 

egi + {H& + H^S* - -^[aVT 7 ?],* + l r Vi = 0. (55) 
Note that for the narrow jet: 

M(r 2 - a 2 ) 

It is also easy to prove that 

{Hv + H^S* = -84,0+4. (57) 
a 

Energy conservation equation in stationary the case may be written 
simply as 

V(a 2 S) = -a 2 a lk r k . (58) 

If we consider that 

(Tiki* = --T^ti « - -T^, r , (59) 
then we can rewrite (58) in the following form: 

V(« 2 5) = aT;/3+. (60) 

It must be noted that our equations contain as "special-relativistic" 
effects (Lorentz factors and all that) as, also, purely gravitational effects 
related to the curving of absolute space (a^l) and the "frame dragging" 
(Pj^O). In the present paper, in purpose to simplify the consideration, we 
shall assume that in the innermost region of the jet, where gravitational 
effects are perceptible, matter moves with non-relativistic velocity (i.e., 
v <Cl and T«l). Thus we shall deal with the general-relativistic but slow 
("non-relativistic" in the sense of special relativity) jet. Hereafter, we shall 
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need to integrate some of our equations over r = const surfaces. According 
to the general theory element of such surface is: 

d 2 x = yJg^geed^dO = a^d(j)d9^(r 2 + a 2 )9d(f)d6. (61) 

First of all, let us integrate continuity equation: 

d r [a^Tmnv r } + de^y^fTmnv 9 =0. (62) 

Remembering that T«l and integrating (62) by d 2 x we get the following 
equation: 

M=2tt f ay^fmnv r d6 = 7r(r 2 + a 2 )9 2 mnv r = const. (63) 

Note that deriving (63) we made the following assumption about the 
angular dependence of n(r,6): 

n(r,0)=n(r)/(e), (64) 

where /(e) > is the dimensionless number-density "profile- function" 
normalized in such a way as to get: 

jj{e)de=l. (65) 

Similarly, under same assumptions, we can integrate toroidal component 
of momentum conservation equation (53). We assume that the jet matter 
may be treated as the medium with ultrarelativistic temperature 

P = nKT, (66) 

e = ran + (7 - iy l P. (67) 
where 7 = 5/3. In this case we have: 

P + £ = mn[l + — — J. (68) 

Taking into account (68), after integration of (53) we get: 
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where L is some constant of integration and £ is defined as: 

C=J Q ef(e)de. 



(70) 



Integration of the energy conservation equation also leads to the another 
algebraic equation of the following form: 



aM 
"2tT 



2/i 

where u;= — (3^ and a) is the another constant of integration. 
If we assume that 5KT/2/j,<^1, then (69) can be rewritten as: 



n = 



ML(r 2 + a 2 )9 2 u; + L 



2\n2, 



2na 



V 



a$ 2 o 



(72) 



For the physical jet solution, the numerator and denominator of (72) 
must vanish simultaneously at a distance r = r A , termed the Alfven point 
of the flow. At this point: 



a 2 A & 



T 



8M£(r 2 A + a?)6\ 



L = 



Mm 



r\ + a 2 )6 A (n - u A ) = 



a A 



'0 — ^A, 



16nV r 



(73) 



(74) 



1-kcla ±un v A 

If we neglect similarly term containing temperature in (71), then it may 
be soled together with (72) for the angular velocity. In such a way we get: 



Q = u + 



a[a+{27rL/M){tt -u)} 



£{r 2 + a 2 )6 2 {n -u) 
Boundary condition Q(ro) = Qq implies that 



(75) 



2ttL 



£(ri + a 2 )9 2 (n - u; o y 



(Q -u) = -ai + 



M a 
and (75) may be rewritten in the following form: 



Q = u + 



a 



[r 2 + a 2 ){n -to) 



a-ap (rl + a 2 )(n -u Q ) 2 ( 6q 
HO 2 " n () U 



(76) 
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Knowing explicit analytical expression for Q(r) we can calculate all other 
physical quantities connected with it. In particular, we can get for toroidal 
magnetic field: 



and then, we can find the expression for radial velocity v r all other jet 
physical variables. 

It must be emphasized that solutions contain unknown "jet profile" func- 
tion 9{r). It may be calculated from the radial component of momentum 
conservation equation, since it leads to the first order differential equation 
for 9(r). 

4 Conclusion 

In this letter we demonstrated the general scheme for the construction of 
the general-relativistic model of the magnetically driven jet. The method 
is based on the usage of the 3+1 MHD formalism. It is shown that the 
critical points of the flow and the explicit radial behavior of the physical 
variables may be derived. All jet characteristics may be expressed through 
one quantity: the jet "profile function" 9{r). The latter quantity may 
be modelled in some way (i.e., by adopting the simplest constant open 
angle jet approximation, or by using of some phenomenological jet profile 
form). Alternative, and more self-consistent, approach should imply the 
solution of the complex first order ordinary differential equation for the 
6(r) function. 

However, the full examination of the problem is beyond the scope of 
this work. The results of this letter are of the preliminary character. Full 
consideration of the problem is in preparation and will be published else- 
where. 



B 4> 



8MA 1 \(r 2 + a 2 )9 2 (n-u) 



$ ) ct{r 2 + a 2 )6 2 _ a 
a 2 e\(r 2 ^a 2 )el(n Q -ujQ) 




(77) 
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